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RESEARCH MEMORANDUM 


SOME EXAMPLES OF THE APPLICATIONS OF THE TRANSONIC 
AND SUPERSONIC AREA RULES TO THE 
PREDICTION OF WAVE DRAG 
By Robert L. Nelson and Clement J. Welsh 

SUMMARY 


The experimental wave drags of bodies and wing-body combinations 
over a wide range of Mach numbers are compared with the computed drags 
utilizing a 2^-term Fourier series application of the supersonic area 
rule and with the results of equivalent -body tests . 

The results indicate that the equivalent-body technique provides a 
good method for predicting the wave drag of certain wing-body combina- 
tions at and below a Mach number of 1. At Mach numbers greater than 1, 
the equivalent-body wave drags can be misleading. The wave drags com- 
puted using the supersonic area rule are shown to be in best agreement 
with the experimental results for conf igurations employing the thinnest 
wings. The wave drags for the bodies of revolution presented in this 
report are predicted to a greater degree of accuracy by using the frontal 
projections of oblique areas than by using normal areas. A rapid method 
of computing wing area distributions and area-distribution slopes is given 
in an appendix. 


INTRODUCTION 


The area rule, first advanced by Whitcomb in reference 1, has con- 
siderably altered the methods for predicting wave drag of wing-body com- 
binations. Studies leading to the discovery of the area rule showed that 
interference drag between wing and body components could be very large. 
Therefore, estimation of drag by component buildup without somehow evalu- 
ating the interference drag could give misleading answers. However, In 
consequence of the transonic area rule, a valuable tool was made avail- 
able to the designers in assessing the transonic drag. This was the 
equivalent -body concept, which states that at transonic speeds the pres- 
sure drag of the airplane Is the same as that for a body of revolution 
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having the same longitudinal distribution of cross-sectional area. As a 
result , the drag of the conf iguration is obtained by either estimating 
or experimentally determining the equivalent -body drag. Experimental 
checks for airplane configurations presented in reference 2 generally 
support this concept in the transonic speed range. 

The supersonic area rule, given by Jones in reference 3, provided a 
powerful method for calculating the wave drag at supersonic speeds. In 
references 4 and 5 the mechanics of the drag calculations were discussed 
together with a number of comparisons of calculated and experimental 
drags generally at low supersonic speeds. Jones pointed out in refer- 
ence 3 that the method could be expected to give good results for thin 
wings mounted on vertically symmetrical bodies. Later, Lomax In refer- 
ence 6 gave the complete linearized theory expressions for the drag. 

The added terms In Lomax's result represented the limitation pointed out 
by Jones. — - 

The purpose of the present paper is to provide a better feel for 
the range of applicability of both the transonic and supersonic area 
rules. For the transonic area rule, this is done by making additional 
comparisons between equivalent -body and wing-body experiments. For the 
supersonic area rule, comparisons are made of calculated and experimental 
results for both body and wing-body combinations over a wider range of Mach 
numbers than heretofore made. The supersonic-area -rule calculations were 
made by using a 24-term Fourier series expression for the slope of the area 
distribution. 


SYMBOLS 


A frontal projection of the area cut by a Mach plane or wing 

aspect ratio 


a 


n 


C D 



drag coefficient, - 2 - 
’ qS 


pressure coefficient 


wing local chord 

root chord of particular pointed wing tip 


wing root chord 
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wing tip chord 


drag 


maximum tody diameter, 2r m 
resultant pressure force 
wing-thickness-distrfbution function. 


wing-area-distribution function 


wing-area-distribution slope function 


/ b cos e\ / 

mil + 4 for the left-wing panel; m l 

\ tan A / \ 

the right-wing panel 


p cos 0 
tan A 


length of configuration 

total length of area distribution 

body fineness ratio 


Mach number 


tan A 


integer 

dynamic pressure 
local body radius 
maximum body radius 


reference area 


body frontal area 


wing exposed area 


wing total plan-form area 
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B 0 

t/c 

x,y,z 

x 0 


max 


wing semi span 

semispan of particular pointed-tip wing 
wing thickness ratio, 

Cartesian coordinates 

point of intersection of Mach plane with the x-axis 

x-coordinate measured from wing leading edge 
dummy variable 

local maximum wing ordinate 


= /m2" 




0 = 


cos 


A 

A 


- 1 ' 2 £ - *) 

wing leading-edge sweepback angle 
wing taper ratio, — 


Mach angle, sin -1 i 


v = 


x' 


value of v at root of particular pointed-tip wing 


value of v at root of actual wing 


angle between z-axis and line of intersection of Mach plane 
with the y, z plane 


REVIEW OF THE BASIC THEORY 


From reference 6, the equation for the wave drag of~any system of 
bodies or wings and bodies can be written as : 
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dx 



dxi 


d 2 A(x,9) 
dx 2 


J3_ df(x,9) 
2 q dx 


d 2 A(xi,9) 



JL ^(x^e) 

2 q dxq 


log(x - Xp) 


( 1 ) 


The equation is subject to the usual limitations of the linearized theory. 


Before discussing the terms in the drag equation, it is well to review 
the definition of Mach planes. The physical significance of equation (l) 
is understood if the configuration is cut by Mach planes. Mach planes are 
easily visualized by considering a Mach cone originating at a point on the 
x-axis which is alined with the remote relative wind. A Mach plane is 
simply a plane tangent to the Mach cone and at an angle of roll, 9 about 
the x-axis measured from the y-axis. By moving the vertex of the Mach 
cone along the x-axis, a series of parallel Mach planes will cut the con- 
figuration for a fixed roll angle 0 . 


In the drag equation the term A(x,9) represents the frontal pro- 
jection of the oblique area cut by a particular Mach plane, whereas f(x,0) 
represents the net force normal to the stream direction on this section in 
the 9 direction. These relationships are illustrated in figure 1 for 
angle of roll 9 of the Mach plane of 0° and 90 °. By neglecting the 


term 


the equation reduces to the supersonic-area-rule formula 


3 df(x,9) 

2 q dx * 

given by Jones in reference Evaluation of f(x,9) requires the pres- 
sure distribution on the configuration which when integrated over the 
configuration gives the drag directly. As a result, large values of 


impose a limitation on the supersonic area rule, even within 
the framework of the linearized theory. 


It Is not the purpose of the present paper to evaluate the drag of 
configurations Including the effect of the pressure term but to evaluate 
the drag of configurations using the supersonic-area-rule formula of 
Jones. The influence of the pressure term was evaluated for one simple 
case. 


Equation (l) can be written In coefficient form as 

°i> = srj 0 c »< e > ae 


( 2 ) 


IBSMTI^L? 
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where 


C d (6) =-^ 


2itS - 


0 0 0 


H 4!a ^Vd^A _ £_ df_\iog(x - xq) 

\dx^ ^ 2q dxj 


dx dx 1 (5) 


The quantity Cj)(0) is most readily determined by solving the inte- 
gral for Cj-)(9) through a Fourier sine series expression_f or ^ fol- 
lowing the method of reference 7 if 


(j> = eos“^ (2 ~ — 1 


a n = | 


’0 


'i - | n 0 d0 

then C-Q ( 0 ) can be written as 

c d( 0 ) = fe Zna n 2 

For the computations of this paper, only 2k terms were used in the 
Fourier sine series expression for Thu6, 


dx 


n=24 


c d( 6) y , *** 


n=l 


BODY DRAG RESULTS 


For bodies of revolution, the calculation of the drag is simplified 
to some extent because the area distributions are identical for all roll 
angles. However, except for high-fineness-ratio bodies, it is not pos- 
sible to assume that the frontal projection of the oblique area cut by 
the Mach plane is the same as the normal area. “Figure 2 shows an example 
of this for two parabolic bodies of revolution having different fineness 
ratios and shapes. The area-distribution curve slopes were calculated 
from the expression 


dA = _2_ f *u (x - xp)dx 

dX p2 ^ x l yp^Cx) - (x - xq ) 2 
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The derivation of this expression is given in appendix A. It has also 
been assumed for the calculations (and all succeeding body calculations) 
that a cylinder can be added at the base of the body without altering 
the drag. If this were not done, the solution would require the flow to 
fill the area behind the base which would exceed the limitations of the 
linearized theory. Figure 2 shows large changes in the peak slope over 
the afterbody of the fineness-ratio-6. (A configuration; these changes 
would lead to a significant drag variation with Mach number. 


The evaluation of the slope of the oblique area distributions is 
extremely difficult except for simple bodies. There naturally arises 
the question as to whether this is worth while if the pressure term Is 
ignored . 

As derived in appendix B, the local force acting on the oblique 
area of a body of revolution is 


f _ 



f dCp 

The only assumption made in derivation of is that is 

constant over the oblique area. This is a reasonable assumption except 
for bodies having discontinuities, and high local slopes. Then, 


dA _ £ f = dA + dC P 
dx 2 q dx 2 dx 


Thus, the error in the drag introduced by ignoring the pressure term is 

dC p 

dependent on the pressure gradient . 

It would be expected that the drag for a conical nose with an 
attached shock wave over which the pressure is constant at zero angle of 
attack would he least affected by the pressure term. (The pressure term 
takes on a value only near the juncture with the cylinder; however, the 
pressure term was not evaluated in this region.) Figure j presents a 
comparison of the drag of various cones calculated with the supersonic- 
area-rule formula with the exact theory drag of reference 8. The lowest 
Mach number of the comparison corresponds to the lowest Mach number for 
entirely supersonic flow on the cone as calculated with the exact theory. 
The highest Mach number of the comparison was arbitrarily taken as that 
at which the slope of the Mach line equaled one-half the slope of the 
nose. The agreement between the two theories is remarkable, within 
5 percent except for a few points . 
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A better comparison may be made by plotting c D (z /r m ) against 
P r m Q. uan 'ti'ty which defines the frontal projection of the oblique 

area distribution. This has been done in figure 4 to give drag in dimen- 
sionless or collapsed form. The drags from exact theory (5° half angle 
cone was chosen as representative) , slender body theory (ref. 9)> and 
supersonic-area -rule theory are shown. The comparison shows the great 
improvement of the area-rule theory over the slender-body theory at 
values of greater than 0.2, and the good agreement of the area- 

rule result with the exact theory to ^r m / , Z of about 0.7* At higher 
values of (3r m /Z the area-rule theory is in error, possibly first because 
the pressure term is neglected but finally, near )3r m /Z = 1, because the 
assumptions of the linearized theory are violated. At (3r m /Z = 1, the 
Mach line lies on the cone surface, which corresponds to the realm of 
hypersonic flows. (See ref. 10.) 

For a body with curvature, for example, a nose of parabolic profile, 
the pressure over the nose is variable, and the influence of the pressure 
term may be significant. Figure 5 presents the drag for noses of parabolic 
profile in collapsed form. Here the supersonic-area-rule theory is an 
improvement over slender-body theory but in only partial agreement with 
the more exact second-order theory of reference 11. Inclusion of the 
pressure term, evaluated by using second-order pressure ~distribut ions, 
however, does give agreement with some of the second-order-theory results. 
Since the second-order-theory drags do not collapse into one curve, agree- 
ment should be. expected only with those points for which the pressure dis- 
tribution used in evaluating the pressure term apply. However, this was 
not the case. For example, the pressure distribution used for the pres- 
sure term calculation at (3r m /Z = 0.3 corresponds to the flagged symbol. 

For the parabolic noses, both the area-rule theory and the area-rule 
theory plus the pressure correction cannot be expected to apply near and 
above (Zr m/ /Z = 0.5, where the slope of the Mach line equals the slope of 

the nose tip. 

Figure 6 presents a comparison of the pressure drag from supersonic- 
area-rule theory with experiment and slender-body theory for a family of 
parabolic bodies of revolution. The experimental drags were taken from 
references 12 and 13 . In determining the experimental pressure drags, 
the friction drag was assumed turbulent and evaluated by using the sub- 
sonic drag level and the results of reference" 14 for the effects of Mach 
number and Reynolds number; the fin pressure drags were assumed identical 
and taken from reference 13 ; and the base drags were small and were sub- 
tracted when available. The slender-body-theory drags were calculated 
using the curves of reference 9« 

mtiS& 
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As would be expected, the comparisons show the increasing ability of 
both the area-rule theory and slender-body theory to predict the drag as 
the body fineness ratio is increased. In most cases, the area-rule theory 
offers a significant improvement over slender-body theory. The area-rule 
theory and slender-body theory are in agreement near M = 1, since at this 
Mach number the supersonic-area-rule theory reduces to slender-body theory. 


From these nose and complete-body comparisons that have been made, 
the following conclusion can be drawn. The area-rule drag of bodies can 
be predicted to a greater degree of accuracy by using the frontal projec- 
tion of oblique areas at a given Mach number than by using normal areas, 
if, at the Mach number under consideration, the limitations of the line- 
arized theory are not exceeded. This is illustrated by the comparison 
between the drag at a given Mach number and the drag near M = 1 especially 
for the low-fineness -ratio bodies. It is not to be inferred from the above 
statement that the supersonic-area-rule method is recommended for evalu- 
ating the drags of bodies of revolution. However, when the drags of wing- 
body combinations for which the body area distribution is needed are deter- 
mined, the oblique area distribution should be used if the body is of low 
fineness ratio or has low-fineness -ratio components. 


CALCUIATION OF WING-BODY DRAG 


The difficulty in computing the wave drag of wing-body configurations 
can be considerably reduced if the configuration meets the following con- 
ditions: first, the body is of sufficiently high fineness ratio so that 

the change in body-area distribution with Mach number is small, and sec- 
ond, the wing is thin. These conditions imply also that the pressure 
term Is negligible. Some feel for the body fineness ratios necessary for 
the above condition to be met can be obtained from the preceding section 
on bodies of revolution. The assumption of a thin wing allows the Mach 
plane intersecting the wing obliquely to be replaced by a plane perpen- 
dicular to the wing chord plane intersecting the wing plane along the 
same line as the Mach plane. Note that, at zero roll angle, the Mach 
plane is normal to the wing chord plane but Is not normal to the wing 
chord plane at any other roll angle for a Mach number other than M = 1. 
Also the angle between the Mach plane and the normal to the wing chord 

plane is greatest and equal to tan - -^ at a roll angle of 90 ° • 

Appendix C presents a sinqple analytical method for evaluating wing 
area distributions and area-distribution-curve slopes. The curves neces- 
sary for evaluating these quantities (figs. 1 6 and 17 ) are applicable only 
to 65 A series airfoils, but similar curves can be made up for other air- 
foil sections. 
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In order to get an Idea of the applicability of the thin-wing assump- 
tion, a calculation has been made of the true area-distribution-curve- 
slope variation for 60° delta wing having an NACA 65 AO 66 airfoil section 
for a roll angle of 90° and a Mach number of 1.4l4. In order to simplify 
the calculation, the wing was approximated by a sufficient number of 
linear-slope elements to define the airfoil section adequately. With 
this approximation the Mach plane intersection with the wing surface was 
made up of straight lines. The expression for the frontal projection of 
the oblique area was then easily evaluated and differentiated to obtain 
the slope. The results of the calculation are presented in figure 7* 
Although the slopes for the upper and lower half wings are significantly 
different, the total slope agrees almost exactly with the slope obtained 
by using the thin-wlng assumption. On the basis of this result, it is 
felt that the thin-wlng solution should be adequate for wings of present- 
day interest . 


For the wing-body combinations of this paper, an additional simpli- 
fication was allowed in the supersonic-area-rule wave-drag calculations. 
Since the tail fins mounted on the models were thin and relatively small 
(see ref. l6), their drags were subtracted as tares. Then, since the 
bodies for all cases were of high fineness ratio (and identical), the 
body-area-distribution-curve slopes were considered independent of Mach 
number, and the changes in the area-distribution-slope curves with Mach 
number and roll angle were due entirely to the wings. As derived in 
appendix C, the area distribution for a given wing (m fixed) is depen- 


dent only on the value of 


0 cos 0 


Thus, the area-dlstribution-slope 


tan A 

curves for the wing-body configuration are dependent only on the value 
P cos 0 


of 

tan A 

configuration. 


Then, from equation (2) and because of the symmetry of the 


p r*/ 2 

c d (m) = f / c D (e)d0 


In order to obtain the wave drag of the configuration, a plot of Cp 
against 0 is required. This can be computed if a plot of Cp against 

' "t.nr^A 6 ' is £ iven > since the angle 0 is known for fixed values of — ^ 
0 cos 0 


and 
0=0 


tan A 
and £. 


tan A 

The configuration drag is simply the average drag between 


For the wing-body calculations of this paper the bodies were identi- 
cal. The body-area-distribution-slope curves are shown in figure 2(b). 

The curve for M = 1.4l4 was chosen as representative for the Mach number 
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range of interest. The wing area distributions and area-distribution- 
curve slopes were obtained by using methods similar to 1 that given in 

(3 0 q g 0 

appendix C. In addition, a limiting value of - — - — = 0.8 was set for 

"Gan 

configurations having blunt leading -edge airfoils. ^Above - = 1, 

the Mach line lies behind the wing leading edge, and the linear theory is 
no longer valid for blunt airfoilsj) 


An example of the wave drag calculation for the most extreme configu- 
ration investigated (60° delta wing, NACA 65 AOO 6 airfoil) is presented in 

figures 8 , 9, and 10. Figure 8 shows nondimensional plots of ^ against 

6 cos 0 

d for various values of -8 . Figure 9 shows the effect on C-n of 

tan A 

the number of terms in the series solution. Except at = 0 and 0.8, 

tan A 

convergence was apparently obtained within 2k terms. Figure 10 shows the 

3 cos 0 

variation of the area distribution drag with *-r 7 — , the variation of 

tan A 

area distribution drag with roll angle, and the variation of the config- 
uration drag with . ^ -. 

0 tan A 


WING-BODY DRAG COMPARISONS 


Figure 11 presents some wave-drag comparisons for wing-body combina- 
tions. The experimental wing-body results were taken from references l 6 
to 19. The wing-body wave drags were obtained in the following manner. 

The friction drags were assumed to be turbulent and were estimated by 
using the results of reference 14 . Base drags and fin pressure drags were 
subtracted using the results of reference l 6 . The equivalent-body drags 
for a Mach number of 1 were obtained experimentally by using the helium- 
gun technique described in reference 2. Tnese models had four scaled 
tail fins. The friction drag was assumed to be the subsonic drag level 
corrected at higher speeds for Reynolds number and Mach number by using 
the results of reference li)-. Base-drag rise and fin-drag rise were not 
evaluated for the equivalent-body models. These quantities, however, 
should be small in the Mach number range where comparison is valid. The 
supersonic-area-rule-theory drags were evaluated by using the method of 
the preceding section. No attempt was made to evaluate the drag with 
the pressure term included. The drag coefficients presented in figure 11 
are based on total wing area. 

The inability of the supersonic-area-rule theory to predict the drag 
near M = 1 is evident for nearly all cases. However, the agreement at 
the higher Mach numbers between the theoretical drags and the experimental 
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wing-body drags is excellent and within the accuracy of evaluating the 
experimental wave drag, except for three configurations. Two of these 
configurations (figs. 11(c) and 11(g)) had 6-percent -thick wings which 
were the thickest wings investigated. The third configuration (fig. 11(h)) 

had a 4^- percent-thick airfoil but with fairly steep wedge components f 

For these conf igurations, a significant effect of the neglected pressure 
term may be possible. As a result , the drags calculated for configura- 
tions having wings of these thicknesses and sections should be viewed 
with caution. 


The comparisons in figure 11 show that the equivalent -body drags 
give a good approximation to the experimental wing-body drags up to a 
Mach number of 1, except for the two configurations having 6-percent- 
thick wings (figs. ll(c) and 11(g)). This result is in agreement with 
reference 2 which shows the validity of the transonic -area rule decreases 
with increasing wing-thickness ratio. At Mach numbers above 1, the agree- 
ment is variable but tends to be consistent with the flatness of the cor- 
responding theoretical curve. That is, as the theoretical drag variation 
with Mach number becomes smaller, the equivalent body gives a better 
approximation of the supersonic -drag level. _ This would be expected, since 
a flat theoretical curve indicates that the variation in area-distribution 
drag with roll angle or Mach number is small. Then the drag for the Mach 
number 1 or roll angle 90 ° area distribution (corresponding to the equiv- 
alent body) is representative of the configuration drag. 

Figure 12 shows the comparison between experimental-configuration 
drag and equivalent-body drag for an airplane configuration. The com- 
parison shows an extreme example, compared with the relatively good results 
of reference 2, of the inability of the equivalent body to predict the 
supersonic -drag level. The equivalent -body drag is approximately 40 per- 
cent low in spite of the low aspect ratio of the configuration. Apparently, 
the configuration tail surfaces cause the area distribution to change mark- 
edly at low supersonic speeds. Below M = 1, the equivalent body gives a 
fair representation of the configuration drag. The drag of the configu- 
ration minus the tail surfaces could probably be calculated to the degree 
of accuracy shown in figure 11. The influence of the tail surfaces, how- 
ever, may be difficult to evaluate. If the horizontal tail supports a 
load when the configuration Is at zero lift, the Influence of the pres- 
sure term may be significant. Although no supersonic-area-rule drag cal- 
culations were made for this airplane, reference 20 indicates that gen- 
erally good predictions of complete -airplane drag can be made. 1 


Subsequent to the preparation of this paper, NACA RM A 56107 has 
been prepared at the Ames Laboratory and presents supersonic -area-rule » 

calculations for a configuration similar to the one shown in figure 12 
but with small differences in area distribution in addition to the 
absence of a canopy. v 
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EVALUATION OF COMPONENT AND INTERFERENCE DRAGS 


As was shown In the preceding section, the supersonic area rule can 
be a useful tool in evaluating the supersonic drag of a wing-body con- 
figuration. In order to assess the efficiency of the combination as a 
whole, however, the effects of the combination on the component drags 
and the interference drag between the components must be known. The 
supersonic area rule provides a valuable method for evaluating these 
effects . 

The supersonic-area-rule equation can, of course, be used in evalu- 
ating the drags of individual wing and body components. This was done 
for a number of bodies in a previous section of this paper. The same can 
be done for isolated wings . An example of this is shown in figure 13 
where the drag of delta wings having 65 A series sections is plotted in 
collapsed form. The area-rule result is compared with a result obtained 
by the method of Beane (ref. 21). The two methods are just two forms of 
the same linearized wing theory. The agreement between the two methods 
is good. 

An example of the effect of the wing-body configuration on wing drag 
Is presented in figure 14. The calculation is for the configuration having 
the closest agreement between the theoretical and experimental drags 
(fig. 11(b)). In this figure, the drag of the exposed-wing panels based 
on total and exposed wing areas Is compared with the isolated wing drag. 
Separation of the wing panels gives approximately a 10-percent reduction 
in wing drag coefficient at Mach numbers above 1.3. As the Mach number 
approaches 1, this favorable effect disappears . This would be expected, 
for at M = 1 the area distributions of the exposed wing panels and the 
isolated wing would be identical if the body were cylindrical. 

Figure 15 shows an evaluation of the interference drag for the same 
configuration. The sum of the calculated body and wing drags is compared 
with the calculated configuration drag. The curves show a favorable Inter- 
ference effect at Mach numbers below M = 1.3. At Mach numbers above 1.3, 
interference drag is, for all practical purposes, zero. Thus for this 
configuration, at Mach numbers greater than 1.3, the only beneficial effect 
of combining the wing with the body comes from the separation of the wing 
panels . 


CONCLUSIONS 


An investigation has been made of abilities of the equivalent-body 
technique and a 2A-term Fourier series application of the supersonic - 
area-rule method to predict wave drag at transonic and supersonic speeds. 
From the theoretical and experimental comparisons made, the following 
conclusions can be drawn: 





Ik 


NACA RM L 56 DII 


1. The area-mile drag of the hodies of revolution presented in this 
report are predicted to a greater degree of accuracy by using the frontal 
projection of oblique areas at a given Mach number than by using normal 
areas . 

2. The supersonic wave drag of slender-ving — body conf igurations can 
be predicted with the supersonic-area-rule formula. For the wing -body 
configurations investigated, the best agreement was obtained for the con- 
figurations employing the thinnest wings. 

3. The equivalent body technique provides a good method for predicting 
the wave drag of certain wing -body combinations at and below Mach number 1. 
At Mach numbers above 1, the equivalent body wave drags can be misleading. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 6 , 1956. 
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APPENDIX A 


AREA DISTRIBUTION SLOPE FOR BODIES OF REVOLUTION 
CUT BY OBLIQUE MACH PLANES 


The area distributions are identical for all roll angles. For 
simplicity a roll angle of 90° will be used in the derivation. 



The frontal projection of the oblique area cut by the Mach plane (see 
sketch (a) ) is given by the equation: 


A 



dz 


From the equation for the Mach plane. 


z is related to x by 


z = 

P 

and 

dz = — dx 

P 


The equation relating y and x is given by 




(x - xq) £ 


y 


^ - (x - xo) 2 


(Al) 
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Then, 


k = \ Mt&to - (x - xq) 2 dx 
P J x z v 

Differentiating the expression for A gives _ 

dA = 2_ f x u ( x ~ x 0) dx 

dX ° ' P 2 4 - (* -xo ) 2 

where x^ and x u are the roots of x = xq - pR(x) and x = Xq 
respectively. 


(A2) 

(A3) 

pR(x), 
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APPENDIX B 

THE NET FORCE ACTING ON THE OBLIQUE AREA OF A BODY OF 
REVOLUTION AT ZERO ANGIE OF ATTACK 


For a body of revolution at zero angle of attack, the net force is 
independent of roll angle. The derivation will be made for a roll angle 



The net force in the 0 direction (the z-direction for a roll angle of 
90°) (see sketch (b)) can be written as 


Since 


dS 



+ dz 2 



and sin a 


V^dy^ + dz 2 


f 

q. 



The pressure coefficient at zero lift is a function of x only. 
The equation for y in terms of x is given by equation (l) of appen- 
dix A as follows: 


y = l^/pSl^x) - (x - xq) 2 
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and. 


dy = 


p2R ( x )ix - ( x - x o) 


dx 


P V P I 2 R 2 (x) - (x - xq ) 2 


Then, the net force can be written as 


f 

q. 


H 


y ° Cp dy\ 


+ 12 

I upper surface \ J 0 




J lower 


surface 


f _ 2 PO 

<1 '"PJx. 


c p[ p%(i, )j| - (x - *o>] ^ £ r*o c pr ?2R(x) s ~ ( * - *q>] 

u Jp 2 R 2 (x) - (x - xq ) 2 


dx + f 


ax 


x l v/p 2 R 2 (x) - (x - Xq ) 2 


f 2 / x u 

q. P 


Cp[p 2 R(x)|| - (x - Xq)] 
Xz y/p 2 R 2 (x) - (x - Xq ) 2 


dx 


Integrating the expression by parts gives 


I " ' I ,/P S? ^ 2r2(x) ' (x ' x ° )£ 
X 1 


If ^£e. is essentially constant between the limits of integrations, 
dx 


h ^ x i 

Then, from equation for the frontal projection of the oblique area 
in appendix A (eq. (A2) ) 

I = -PA ^ 

q dx 
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APPENDIX C 


METHOD FOR DETERMINING WING-AREA DISTRIBUTION AND 
AREA-DISTRIBUTION-CURVE SLOPE 


This method assumes that the wing is thin and that the oblique Mach 
plane can be replaced by a plane perpendicular to the wing chord plane. 
The method also assumes the wing has straight leading and trailing edges 
and constant thickness ratio. 

The method is developed first for pointed-tip wings. Then, correc- 
tions are made for curved-wing — body junctures arjd finite wing tips. In 
addition, the right- and left-hand wing panels are considered separately. 


Pointed-tip wings .- Consider the right-wing panel shown in the fol- 
lowing sketch: 



(c) 

^^p^IDEiro^ 
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The frontal projection of the area of one wing panel cut by the Mach 
plane is given by 


A = 


2 z 


dy 


and — can be written as 
c 


z fmax _z_ = 1 t , v 

c c z max 2 c ^ ^ 

Then, for — constant, 

A = s o c o I / f(v) 

The value of r) and v are related by the intersection line of the Mach 
plane and the wing chord plane for the right-wing panel given by the 
equation 

x' = Xq + y(p cos G - tan A) 
and for the left-wing panel by the equation 


Then, 


With 


x' 


= Xq + y(p 


cos 6 + tan A) 



+ 


^0 

c 0 


(p cos 9 ± tan A 


h 


c 

c 0 


= 1 - 


n 


V 



+ 



cos 0 ± tan A)t) 


and 


T) = 


v 0 - v 


(tan A i p cos 0)— - v 
c 0 


Let 


m 



tan A 
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In terms of tapered wing geometry, m is given by 


m = A 2- , ±. JN tan A 
A 1 - A 


Then, 




v Q - v 


m 


(l ± P 003 Q ) - v 
V tan A / 


Let 


K 


= m 



3 cos 9 
tan A 


for the left-wing panel and 


K 



(3 cos 9 
tan A 


for the right-wing panel. Then, 



£2. K ~ v 0 

dv (K - v ) 2 


c , - K - v 0 

£7 ~ 1 n " 


(k - Vq V 

dTj = -1 dv 

(K - V ) 5 

The equation for the area can then be written as 


c 

c 0 


a ■ -0% 1 (k - v o) 2 av . 


dower < K - v > 3 
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The slope of the area-distribution curve is obtained by differentiating 
the expression for A 


dA 

dVr 


8 0 c 0 I 


l(K - V 0 ) 


‘'upper f ( v ) 

(K - v)' 


dv + 


lower 


(K - vp) 2 
~ v lower) 


^■ v lower f. / v \ ^ 2^ 

’ 1 flower ) “ 


dv 


upper 


dv 


0 


K - v \ 5 dv n 
upper) 0 


f ( v upper) 


or 


dA 

dv 


'0 


= s Q c 0 | H(K,vq) 


Curves of G(K,vo) and H(K,vo) have been made up for a 65 A series 
airfoil and are given in figures 1 6 and IT for values of K from 0 


to 2.b and v 0 from 0 to 1. In evaluating 


f(v) 


dv, f(v) was 


(K - v) 3 

assumed to vary linearly between airfoil ordinate stations* Figure 18 
gives a plot of f(v) for this assumption. 

For K and v Q greater than 1, G(K, Vq) and H(K, v 0 ) are given 
by the expressions 

G(K,v 0 ) = G(K,1) (K ~ VQ) - 
(K - l) 2 


H(K,v 0 ) = H(K,1) 


(K - vq) 
(K - 1) 


For K and v Q less than 0, G(K,v 0 ) and H(K,vo) are given by 

the expressions 


G(K,v 0 ) = G(K,0) 


K c 


h(k,v 0 ) = H(k,0) (K ~ V °^ 

K 
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Correction for curved-wing — body juncture .- The following sketch 
shows a pointed-tip wing mounted on a curved body. 



(a) 

The areas and slopes will be referred to the actual wing geometry 

(c r , s, and A) . The areas and slopes, however, will be for the exposed 

pointed-wing tip. 

In sketch (d) consider one point of intersection of the wing panel 
with the body. The area of the wing cut by the Mach plane through this 
point is determined only by the product of s 0 c Q of the exposed wing 

through the point and the value of v 0 for the exposed wing at the point 
of intersection. As the point of intersection changes, s Q , cq, and Vq 

change and account for the intersection line. Expressed in terms of the 
actual wing-body characteristics, SqCq is given by 



The quantity £ is related to v n by the expression 

S \J 

v - (1 - A) m I 

v = 

1 - (1 - A)f 

The area of the exposed wing panel cut by the Mach plane can be written 
as 


IDENT^A'j 



2k 


It 
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A = 3-^(1 " (1 ” A )|] G ^ v 0) 

The area is calculated for given value of Vq. The center-line value of 
v is given by 

v r = v b [l - (1 - A)|] + K(1 - A)§ 

The slope is obtained by differentiating the expression for A and 
is given by 


c r s £[l - (1 - A)|] | 

[l - (1 - A)|]H(K,vo) - 2(1 - A)^| G(K,v 0 ) 

1 - A 

1 - (1 - A)| + (K - v 0 )(l - A)ll 


d £ 

if — a = 0 . 

dVQ 

(i-a#(k,v 0 ) 


Correction for finite wing tip .- In order to correct the pointed- 
tip wing panel and slopes for the finite wing tip, the areas and slopes 
outboard of the wing tip are subtracted. 



From sketch (e): 


s o c o = 


c r s 




1 - A 
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Then the areas and slopes are given by 


C r S 


M;ip _ 


1 - A 


Gfcv 0 ) 


^ip c r s r * 


dv 


h(k,vq) 


' r 1 - A 
The center -line value of v is given by 


v r = v q a + K(1 - A) 
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£xecf theory Ah =5°J 
Area- rule theory 
S lender- body theory 


Figure 4.- The drag of cones in collapsed form 


S tender- foody theory 
Area-ru/e theory 
— Area-rule theory + pressure tenr> 
O Second-order theory 


Figure 5>- The drag of parabolic noses in co ll apsed form. 






1/d • 6.04 


l/tt - 8.9/ 


1/d^/tS 


(a) ~ = 0.4, 0.6, and 0.8. 

Figure 6.- Comparison of the drag of parabolic bodies calculated with the 
area rule with experiment and slender-body theory. 










Afatfr p/cme 
jnfersecC/on 

Upper surface- 
lower surface— 


-L«A- o 

s.tp dx 


-Lower naff mng 
-Upper Pa If m'/?q 
O Thin wing 


Figure 7.- Comparison of the area-distribution- curve slopes for a delta 
wing obtained by using the thin- and thick- wing assumptions. MCA 
65AOO6 airfoil section; M = 1 . 4 l 4 ; 6 = 90 °. 
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Figure 8 .- Example of the area-distribution-curve slope for a wing-body 
configuration for various values of p cos 9/tan A. 60° delta wing; 
NACA 65 AOO 6 airfoil section. 
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(a) Area distribution drag. 



(b) Area distribution drag against roll angle. 



(c) Configuration drag. 

Figure 10.- An example of the calculation of the configuration drag from 
the drag of the area distributions at various values of (3 cos 9/tan A. 
6o° delta wing; NACA 65 AOO 6 airfoil section. 
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.O! 


2 Exper/ment 

■ Areo-ntb theory 


=3 


/■O A/ f.2 AS 74 /.S /,6 

/W 


(a) Basic body; model 1 (ref. l6) ; ^ = 0.0305. 

Sw 



(b) Model k (ref. 1 6) ; 6o° delta wing; MCA 65AOO3 airfoil section; 


Sb 

— = 0.0305. 
Sw 



(c) Model 5 (ref. l6); 6o° delta wing; MCA 65AOO6 airfoil section; 

Sb 


Sw 


= 0.0505. 



(d) Model 6 (ref. 19); 6o° delta wing. Thickness ratio varies from 0.03 

Sb 

at root to 0.06 at 0.9 semispan. — = O.O305. 

Sw 

Figure 11.- Comparison of the calculated drag with experiment and 
equivalent- body test results for wing-body combinations. 
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(b) Vq from 0.5 to 1.0 
Figure 17-- Concluded. 

AL* 
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Figure l 8 . - Approximation of 65A series airfoil for the calculation of 

G(K, Vq) and H(K, Vq) . 
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